Phase Diagram for Quantum Hall Bilayers at v = 1 
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We present a phase diagram for a double quantum well bilayer electron gas in the quantum Hall 
regime at total filling factor v = 1, based on exact numerical calculations of the topological Chern 
number matrix and the (inter-layer) superfluid density. We find three phases: a quantized Hall 
state with pseudo-spin superfluidity, a quantized Hall state with pseudo-spin "gauge-glass" order, 
and a decoupled composite Fermi liquid. Comparison with experiments provides a consistent expla- 
nation of the observed quantum Hall plateau, Hall drag plateau and vanishing Hall drag resistance, 
as well as the zero-bias conductance peak effect, and suggests some interesting points to pursue 
experimentally. 
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The coexistence of an incompressible integer quantum 
Hall effect (IQHE) state and interlayer superfluidity has 
been established through a series of experimental and 
theoretical works on bilayer two dimensional electron sys- 
tems at a total electron filling number v = liSiiii. Such 
an IQHE is a consequence of the strong Coulomb in- 
teractions, which lead to a charge gap at v = 1. De- 
noting the layer index as pseudospin "up" and "down" , 
the ground state is a quantum ferromagnet with spon- 
taneous interlayer phase coherence^, which exhibits su- 
perfluidity in the zero layer separation limit, d = 0. On 
the other hand, in the d — > oo limit, the two layers 
are decoupled, each comprising a compressible compos- 
ite fermion liquid at v — 1/2. The nature of the phases 
and transitionsi&SiiSiii between these two limits has at- 
tracted many recent studies. A generalized pseudospin 
description suggests a first order transition^ at which the 
pseudospin order vanishes upon increasing d. The rela- 
tion of the loss of ferromagnetism with the disappear- 
ance of the IQHE remains unclear, and indeed several 
scenarios in which the ferromagnetism and IQHE do not 
vanish simultaneously have also been proposed based on 
the Chern-Simons mean-field theory^. Furthermore, in 
real samples with impurities, a first order transition is 
believed to be impossible on general grounds 1 ^. A con- 
sistent picture for the phase transitions is still absent. 

Fundamentally, the distinct phases of the system at 
T — are characterized by their topological order and/or 
broken symmetries. Experimentally, these are predomi- 
nantly reflected in electrical transport coefficients, partic- 
ularly the Hall conductance. Due to the lack of exact so- 
lutions (except at d = without impurities), the ground 
state phases and transport properties have mainly been 
discussed based on effective theorie a 5 ' 7 ' 9 . Exact numeri- 
cal calculations for these systems have been done in the 
absence of random disorder potentiali2*ii, which cannot 
provide direct information about the transport. In this 
letter, we report the first finite-size exact calculations of 
transport properties for such a system in the presence of 
random impurities, by obtaining the topological Chern 



number matrix of the many-body wavefunction and the 
superfluid density of the ground state. 

The Chern numbe»i2*i4 is a unique integer topologi- 
cal invariant associated with a wavefunction, and can be 
used to distinguish different quantum Hall states. Phys- 
ically, the Chern number equals the boundary condition 
averaged Hall conductance (in units of e 2 /h), so an IQHE 
state is expected to display a fixed non-zero integer Chern 
number independent of disorder configurations. A state 
with a non-quantized Hall conductance instead displays 
a random integer Chern number intrinsically fluctuating 
with different disorder configurations or other external 
parameters. Such states are generally expected to be 
critical or fluid in nature, since a nonzero currenlii 4 * 1 ^ 
necessarily exists in the bulk to destroy the exact quan- 
tization of Hall conductance. Thus the distribution of 
Chern numbers over samples also reveals the extended 
or localized character of the state. 

Numerical calculations of the Chern number have up 
to now only been carried out from single-particle wave- 
functions in non-interacting system o 14 ' 15 . The Chern 
number of a many-body wavefunction is, nevertheless, 
well-defined, albeit difficult to calculate. We have de- 
veloped a new approach to obtain the exact Chern 
number by numerically evaluating the Berry phase of 
the many-body wavefunction upon changing boundary 
phases adiabatically 16 . In the present bilayer system the 
topological Chern numbers form a 2 x 2 matrix related 
to the topological ordering of the systemic, which deter- 
mines the charge and spin (we will refer pseudospin as 
spin) Hall conductances (as well as the Hall drag, which 
is the difference of the two). 

Our main results can be summarized as a numerical 
phase diagram (see Fig. 1) in the d — W plane (W is 
the disorder strength) with three distinguishable phases. 
Phase I is the usual bilayer ferromagnet, embodying co- 
existence of the v = 1 IQHE and inter-layer superflu- 
idity, occurring at small d and W. At relatively strong 
W but small d regime, the superfluid state will first un- 
dergo a phase transition to a "gauge glass" (phase II), in 
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FIG. 1: The solid line is the critical layer separation d c vs. 
W, which separates the IQHE plateau (phase I and II) from 
the composite Fermi liquid state (phase III). The dashed line 
is the critical d s vs. W, which separates the superfluid state 
(phase I) from the gauge glass state (phase II). 

which n s vanishes due to strong phase frustration. Both 
phases I and II have the same Chern number matrix: a 
uniquely quantized charge Chern number C c — 1 and 
a random spin Chern number, which result in nonzero 
Hall drag conductance. In phase I, spin superfluidity 
implies that the full spin resistivity tensor and hence 
p s xy vanishes below the Kosterlitz-Thouless temperature, 
so we expect quantization of the Hall drag resistance, 
Pxy = Pxy ~~ Pxy = h/e 2 with exponentially small ther- 
mally activated corrections at low temperature. Phase II, 
despite having the same Chern number matrix as phase 
I, has a different behavior of the spin resistivity, since the 
average n s vanishes. Conventional theory for the gauge 
glass carried over to this problem suggests^- p xy ~ p xx 
are exponentially small but non-zero at any T > and 
vanish only at T = 0, leading to a weaker quantization 
of the Hall drag resistance. With the increase of d or W, 
phase II undergoes a transition to the compressible Fermi 
liquid (phase III). Phase III can be understood from the 
large d limit of decoupled layers. It is a metallic state, 
characterized by zero drag conductance and a finite di- 
agonal longitudinal conductance. These results provide 
a consistent understanding for the existing experiments 
(see below). 

In the presence of a strong magnetic field, the Hamilto- 
nian, projected onto the lowest Landau level, for a double 
layer two-dimensional electron gas can be written as: 

H = E E e " g2/2 ^(«) eiq ' (R ' _Rf) 

i<j.a,{3 q^O 

+ (i) 

where R" is the coordinate of the i-th electron in layer a 
(a = 1, -1). V a> p(q) = 2ire 2 /(eqL x L y ) x exp(-qd6 a ,-p), 
is the Coulomb potential. VJ" lp (g) is the impurity po- 
tential generated according to the correlation relation < 

Kmp(9)^im P (-V) >= ift"<W<W 5 which corresponds 



to < V& >= W 2 S a , p S(r - r') in real space, 

where W is the strength of the disorder. We set the 
magnetic length I = 1 and the interaction strength 
e 2 /el = 1 for convenience. We impose the generalized 
boundary conditions: T a (L 3 )\$ >= e l9 ?\Q > (J = x,y), 
to the finite size double layer system, each in an L x X Ly 
rectangular cell with an integer number of flux quanta 
N s = L x L v /2tt. T Q (r) is the many-body magnetic trans- 
lation operator. The tunneling term is not considered 
here in order to study the interesting limit where the 
correlation between two layers is purely due to Coulomb 
interaction. We consider up to N e — 12 electrons at fill- 
ing number v — N e /N s — 1, spanning a Hilbert space of 
size Nbasis = 853776. 

Through a unitary transformation "F = 

ezpH£iiiEa(£z? + zjw?)]*. the topological 
Chern numbeniii^ can be calculated as: 

i P <9<F 

^-5 •!*?>- <Sjl»»- (2) 

where 8j has layer index a and (3 with 9j = O",0^. With 
a, j3 = —1,1, C a ^ forms a 2 x 2 matrix. The closed 
path integral is along the phase boundary of the 27r x 27t 
unit cell. If common or opposite boundary phases are 
opposed on the two layers, then one obtains charge and 
spin Chern number, C c,s , which is related to the bound- 
ary phase averaged charge and spin Hall conductances 
a X y = C c,s ^-, respectively. We separate the phase space 
into a mesh of 64-200 squares. By repeatedly calculat- 
ing wavefunctions at all nodes of the mesh using Lanczos 
method, we determine the integer Chern number for the 
many-body state in each disorder configuration. 

We first consider the charge Chern number C c as a 
function of d. At a weak disorder strength W = 0.16 
and d < 1, we find C c = 1 for all disorder samples at 
7V e = 6,8,10,11 and 12 (20 samples for iV e = 12 and 
1000 for N e — 6). Hence the corresponding ground state 
displays the IQHE with total (charge) Hall conductance 
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a xy = Sj-. As we further increase d, a strong fluctuation 
of the Chern number takes place at d ~ 1.1 — 1.3, which 
are caused by level crossings upon tuning the boundary 
phases and disorder. The persistence of the crossing of 
low energy states at particular values of d, W is a sig- 
nature of a quantum phase transition, and is associated 
with the collapse of the mobility gap. Thus the mobility 
of the state in the charge channel is tied to the fluctua- 
tions of C ° 14 ' 15 , motivating us to define p ext — P(C C ^ 1) 
(the probability of finding C c ^ 1) as a characterization 
of the extensiveness of the many-body state. If p ex t ex- 
trapolates to a nonzero value in the thermodynamic limit, 
it represents a fluid phase which can carry current in the 
bulk to spoil the exact quantization of a xy . 

To determine the critical d c for the v = 1 IQHE to 
metal transition, we plot p ex t as a function of d at W = 
0.16, N e = 6, 10 and 11 in Fig. 2. p ext always increases 
rapidly around d ~ 1.1 — 1.3 and saturates to 0.5 at larger 
d. In the same figure, we plot Ap ext /Ad (Ad — 0.05) as 
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FIG. 2: The effective charge mobility p ext at W = 0.16, 
defined as the ratio of the number of samples with C c 7^ 1 
and the total number of samples, as a function of d at N e = 6, 
10 and 11. Ap ext / Ad vs. d is also shown. 

a function of d. Around d = 1.25, there is always a very 
strong peak growing with N e . We determine the critical 
layer separation d c (= 1.25 ± 0.05 for W = 0.16) for the 
IQHE to metal transition by the location of the maximum 
of this peak for each W. A step jump of p ex t is expected 
for such a transition in the thermodynamic limit. 

The critical d c as a function of W determined in this 
manner is indicated in Fig. 1 by the solid line. The d c 
for two aspect ratios (L x /L y = 1 and 1.6) agree well 
with each other, demonstrating the smallness of overall 
finite size effects. At W — > 0, d c saturates to around 
d c = 2.2 ± 0.2 with relatively large error bars due to 
increasing finite-size effects at very weak W (at W — 
0.04, d c = 2.02 ± 0.13). We note that studies based on 
the local ferromagnetic moment tufm for weak disorder 
will tend to underestimate d c , due to a large reduction of 
itifm by low energy states mixing around d w 1.5, which, 
however, does not affect the Chern number. 

More information on the nature of spin sector can be 
obtained by tuning the boundary conditions in two layers 
according to 6t = X = — 0~ l . The energies of the lowest 
two states E g and E\ change significantly with 9 t , with 
variations of the order of the level spacing. At small d 
and W, Eg first increases as a quadratic function of 9 t 
until it is energetically favorable for a vortex (through 
the hole in the torus encircled by the x axis) to enter 
the system, a typical feature of a superfluid state. The 
disorder averaged superfluid density can be calculated as 
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In Fig. 3, n s is plotted as a function of d at W = 0.16 
for N e = 6,8,10 and 12. The overall behavior agrees 
with the generalized mean-field calculation for pure sys- 
tem (our definition of n s is twice that of Ref. @). Nat- 
urally, n s reduces with increasing d, and we define the 
boundary of the superfluid state d s by n s — 0, e.g. 
d s = 1.2 for W = 0.16 shown here. The critical d s 
is also shown in Fig. 1 as the dashed line. In the 
strong W case, d s becomes obviously smaller than d c , 
indicating a superfluid state (phase I) to phase II tran- 
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FIG. 3: The superfluid density n s 
transition point d a = 1.2, n s —* 0. 
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sition inside the IQHE regime. In phase II, the aver- 
aged superfluid density vanishes, due principally to phase 
frustration: E g (6 t ) still depends strongly upon 6 t with 
positive or negative curvatures depending on disorder 
realizations. We postulate hence that the spin sector 
in phase II behaves as a gauge or vortex glassi^, with 
Edwards- Anderson magnetic order and algebraic stiffness 
(\Eg(9 t = tt/2) - E g (6t = 0)|) ~ L"l e l, as proposed for 
v = 1 bilayers in Ref. |5j. 

Further evidence that phase II is not a spin insulator 
is obtained from the spin Chern number C s , determined 
by imposing opposite boundary phases to both layers. 
We find that C s fluctuates (around 1) throughout the 
d — W plane, implying fluidity of the spin sector and a 
non-quantized spin Hall conductance. This rules out a 
spin insulator for phase II, favoring the interpretation as 
a gauge glass. In both phase I and II, the Hall drag con- 
ductance cr^+cr-y 1 ' 1 = \(sJ% v -o s xy ) = {C c -C s )^, is 
also non-quantized and non-zero, which is a consequence 
of the coupling between two layers. At the phase bound- 
ary for the IQHE (d = d c ), we find that the nondiagonal 
Chern number C 1,_1 = C -1,1 = and the drag Hall con- 
ductance drops to zero, indicating that the spin sector is 
also involved in the phase transition at d — d c . At this 
gauge glass to composite Fermi liquid critical point, we 
expect the spin correlations go from (Edwards- Anderson- 
)superfluid to metallic. 

To reveal the charge plateau, we calculate o c xy as a 
function of v. Shown in Fig. 4, at W — 0.16 and d = 
0.8, we observe an exact quantized plateau between v 
111 2 = 0.91 and v = 13/12 = 1.09, with p ext = 
corresponding to a finite mobility gap. The plateau width 
is usually smaller or around Ais — 0.2 depending on W, 
in good agreement with experiments^. In contrast to the 
plateau in charge channel around v = 1, the spin Chern 
number fluctuates with different disorder configurations, 
as shown in the inset of Fig. 4. 

We conclude with some comparison to experiments. 
Both phases I and II exhibit the IQHE in the charge chan- 
nel, i.e. p% y = h/e 2 and p c xx = at zero temperature, and 
we expect activated corrections at T > 0. The spin Chern 
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FIG. 4: The charge Hall conductance a% y (in units of e 2 /h) 
vs. filling number v — N e /N s . In the inset, the random spin 
Chern number C s vs. disorder configuration I for 20 samples 
at v = 1. W = 0.16 and d = 0.8 (in phase I) for both data. 

number is random and fluctuating, which gives rise to 
a nonquantized spin conductance. The exact quantized 
Hall plateau (for the charge) has been observed by several 
experimental groups^^ with a critical d e c xp — 1.8, cor- 
responding to the regime of very weak W ~ 0.07 in our 
phase diagram. In this regime, the width of phase II is 
extremely narrow if non-zero, and the Hall plateau phase 
at d < d c observed experimentally therefore corresponds 
to phase I. As a true two-dimensional superfluid, it is 
expected to exhibit anomalous properties, e.g. divergent 
a xx even at non-zero temperatures below the Kosterlitz- 
Thoulcss temperature, T < Tkt- The experimentally 
observed zero-bias tunneling conductance peak below a 
characteristic temperature^ is a direct reflection of the 
associated off-diagonal order in the spin channel. Fur- 
thermore, the divergent spin conductivity implies that 



the full spin resistivity tensor must vanish, p s xx = p s xy — 
for T < T KT . At T = 0, the IQHE then implies quan- 
tized Hall drag p xy = p% — p s xy = h/e 2 and vanishing 
longitudinal drag, p xx — p c xx — p s xx = 0, in agreement 
with experiments^ -. While the quantization of Hall drag 
is expected to be violated for T > by activated pro- 
cesses contributing to p xx , the theory predicts that the 
spin Hall resistivity, should vanish even at non-zero tem- 
perature for T < Tkt, which would be interesting to ex- 
plore experimentally For d > d c , in phase III, we obtain 
C 1 ' -1 = C -1 ' 1 = numerically, signaling the decoupling 
of two layers. As a consequence, the drag conductance 
and resistance reduce to zero at T = 0, which also agrees 
with the experiment at larger d (d > df* p ). Lastly, for 
the relatively pure sample used in Refl20l phase I is di- 
rectly neighboring phase III which results in p xx ^ only 
along the phase transition line d = d c at very low tem- 
perature, a property again observed experimentally 20 . 
For more disordered samples at intermediate d, phase II 
(gauge glass) intervenes. While we expect this phase ex- 
hibits the same transport coefficients as phase I at T = 0, 
it has no associated Kosterlitz-Thouless transition. This 
implies that p xy ,p xx are generally non-zero for T > 
in this range, and probably that p xx is enhanced at low 
but non-zero temperatures. We leave a more detailed 
investigation of the gauge glass to a future study. 
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